Introduction. — The concepts Scaling and Universality have played an essential role 
in the description of statistical systems |l|]. 

Recently a multisite interaction system on Husimi tree approximation was investigated 
0. First, it was shown, that this approach yields good approximation for the phase 
diagrams, which closely match the exact results obtained on a Kagome lattice 0. Second, 
a multisite antiferromagnetic interaction was studied and interesting connections with the 
area of dynamical systems was made. The qualitative picture of full doubling bifurcations 
^ ' diagram including chaos, period-3 windows, etc., for the magnetization of the base site 

^ I of this system was exhibited, whereas in antiferromagnetic Potts model only one period 

doubling occurred 0]. 

On the other hand, it is well known, that universality of Feigenbaum constants directly 
^ ! applies to period doubling bifurcations sequence 0. 

The aim of our paper is numerical calculation of the Feigenbaum constants for the three- 
site antiferromagnetic interaction (TSAI) Ising spin system and to receive the quantitative 
description of the transition of this statistical physical system from ordering to chaos. 
>■ ■ Husimi tree and recursion relation. — The pure Husimi tree ||^ is characterized the 

^ ! 7-the number of the triangle neighbors. The Oth-generation is a single central triangle. 

Q I The TSAI model in the magnetic field defined by the Hamiltonian 



O^ 



m 



o 



Qi H = -j'^^aiajUk-fi^ai, (1) 

^^ ■ A i 



Oh! 



where cxj takes values ±1, the first sum goes over all triangular faces of the Husimi tree 
and the second over all sites. Besides we denote J3 = (3J^, h = (3h , (3 = 1/kT, where 
5J^| h-external magnetic field, T-temperature of the system and J3 < corresponds to the 

i-^ ' antiferromagnetic case. 

The partition function will be written as 

Cd ; Z = ^explj3^cXiajCXk + h'^cXi 

{a} [ A i 

where the summation goes over all configurations of system. 

The advantage of the Husimi tree introduced is that for the models formulated on it, 
exact recursion relation can be derived. When "cutting apart" the Husimi tree at the 
central triangle it separates into 3 identical branches and each of them contains 7 — 1 
branches. Then the partition function mey be written 

Z=Y.exp ( JsE^'^^?^^'^ + hj:4A [9n{4'^)r\9n{a^^^)r\gn{<rl'Y~\ (3) 
{<^o} I A i J 

where cxo are spins of central triangle, n- number of shells and the equation for one of 
branches can be written: 

{<7,^ao} I A A i ) 



One of branches, in its turn, can be cut on the site of Ith-generation, which is the nearest 
to the central site. Therefor, the expression for gn{<^o) can be rewritten in the form: 



gn{<yo) = exp I JsXl^ocrif^i + h^ai\ [gn-ii<y[^^)] [fi-n-ilo-f^)] • (5) 



From eq.(^ one can easily obtain: 

gni+) = e''-'''gr-li+)9Vii+) + 2e-'^g:zli+)9r-lH + e'^-'^gViHgrAH, 

guH = e-''-'''g:zli+)9r-li+) + 2e'^glZli+)gViH + e-''-''gViH9ViH- 
Let the following variable be introduced: 

X„ = ^y (6) 

Then for Xn we can obtain the following recursion relation: 

Xn = fiXn-l), fix) = ^2^2(7-1) +2z^a;^-l + l' ^^^ 

where z = e^'^^, /i = e^'' and < x„ < 1. The eq.(0) coincides with that obtained by 
Monroe 0, when pair interaction absents. 

For magnetization of the central base site we obtain: 

m-lrr)- ^'^^(+) " ^''dZH _ e'xl - 1 

" ^ °^ " e'^g2{+) + e- V(-) " e^^xl + 1 ^^' 

TSAI system and Feigenbaum constants. — As it is mentioned in Introduction of this 
paper, the TSAI system is the nonlinear dynamical system and the qualitative picture of 
full doubling bifurcations diagrams, chaos etc., for the magnetization of the base site of it 
was existed f^. 

The questions we want to address in this paper are, how calculate the constants of 
Feigenbaum for TSAI system and if calculated values will coincide with the famous uni- 
versal Feigenbaum constants: 

a = 2.500290 ... , 5 = 4.669201.... (9) 

Feigenbaum observed for logistic map (see ref.[^) two kinds of scaling: one that the 
length 2" cycle first appears at a r„ value, which obeys: 

Tn = i^oo — const6~"', n ^ 1, (10) 

where Too the value of r from which the chaotic behavior ensues and the sequence essentially 
never repeats itself. 



The other scahng was a special behavior which occurred near the x* -value for which 
the map is extremal {x* = 1/2 in the logistic map). If one started out at value for x* then 



a = ^^, n>l, (11) 



tin 
dn+1 

where 

dn = fl~\x*)-X*. (12) 

In eq.(|l2D Rn are the values of r (ri < i?i < . . . < i?„ < r„) and 

flix*) = x*. (13) 

Note, that values of Rn and r„ have the same scale and Too = -Roo- Therefore 

Roo = Rn- constS^''. (14) 

Eqs.(0) and ([I^ defines two Feigenbaum constants, which turns out to be "universal". 

Now let us turn to our questions. One can see from eq.(|^), that the role of above 
mentioned parameter r for TSAI system on Husimi tree for each fixed temperature plays 
external magnetic field h. The recursion function of eq.(0) has one maximum at x* = 
1/ ''-^JI. Note, that this x* depend on values of T and h, whereas in case of logistic map 
it is a constant. Further, we numerically solve the eq.(|l3]) and find out the values of 
Hn {Hn is the analog of i?„). Using this values of Hn and eqs. ([TT|) and (|1^ we calculate 



the Feigenbaum constants for TSAI system. All our numerical calculations are done for 
7 = 3, T = 0.3 and J3 = —1, and are listed in table 1. 

For const., presented in eq. (p!4D , which is depend on family of reflection functions 0, 
for TSAI system we obtain the following value: const. = 0.99 . . ., whereas for logistic map 
it is 0.12.... 

Using the values of Hn and corresponding them xi,X2, . . .Xn, 'we can also calculate the 
magnetization for base site (for each cycle of period doubling) of this system by eq.(|D. 
It means that each 2" period doubling have n values of magnetization, which should be 
explained as an arising of a n-sublattice phase such that xi,X2, . . .Xn determine the states 
on each sublattice. 

One can see from table 1, that for real statistical physical system obtained values of a 
and f5 coincide with famous Feigenbaum constants (eq.(^) with high accuracy and thereby 
confirm they universality once more. 

The some numerical values for magnetization of base site of TSAI system are listed in 
table 1. as well. 

It is interesting to note, that if one lets 7 = 2 in eq.(0) rather then 7 = 3, the above 
mentioned situation changes dramatically. 

Let us consider the following system of equations: 

fix) = -1 ^^^^ 

The eq.([15|) determined the point, where the first doubling bifurcation is begun. 



For recursion function when 7 = 2, the eq.(|T5|) will have the form: 

j n'^x^ + zfi{2 — fi)x'^ + (1 — 2;u)x — z = , . 

\/iV-2Vx-(l + 2/x) = ' ^ ' 

which for any T and h have only nonphysical solutions. Therefore, for TSAI system when 
7 = 2 the period doubling bifurcations picture absents. It means that for this statistical 
physical system there is not phase transition of second order when 7 = 2. 

Conclusion. — In this paper we have investigated TSAI Ising spin model by approx- 
imating it with Husimi tree structures and calculate the Feigenbaum constants a and 6. 
The numerical results show, that obtained valued for these constants for real physical 
system coincide with the famous universal Feigenbaum constants with high accuracy. The 
quantitative description of the transition from ordering to chaos is also obtained. Hence 
we see many of the very intensely studied and by now familiar properties of dynamical 
systems theory, which gives possibility to study the statistical physical systems in a new 
context and in a simple manner. In particular, in this paper with using the well known 
technique for dynamical systems, we analitically show, that for TSAI system there is not 
phase transition of second order when 7 = 2. 

We think, that obtained results are interesting and we plan to continue to investigate 
this line of approach for TSAI system and for several other systems. 

On the other hand, the study of chaotic statistical physical system has opened new 
challenges for theories of stochastic processes, especially in the direction of stochasticity 
of vacuum in QCD 0. In this direction the interesting results for Z(Q) gauge model 



on generalized Bethe lattice was obtained [10|. It gives bases to suppose that TSAI 



Ising spin model on Husimi tree approximation can be connected with double plaquette 
representation of the gauge theory. 
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Table 1. 



period 
doubling 


Hn 


a 


5 


magnetization 
m 


^n 


2^ = 2 


0.18354515... 






-0.6782977 
0.09151673 


0.5423560 
0.999999 


22 = 4 


0.28692571 . . . 


4.86428158 . . . 


3.50752342 . . . 






2^ = 8 


0.31861247... 


2.19505287... 


4.32097441 . . . 


-0.8924331 

-0.1373881 

-0.9581313 

0.1182513 

-0.8506408 

-0.2286943 

-0.9640809 

0.1579719 


0.3457500 
0.8200590 
0.2495890 
0.9733620 
0.3886090 
0.7699640 
0.2369180 
0.9999 


2^ = 16 


0.32607381 . . . 


2.76000232 . . . 


4.5870529 . . . 






2^ = 32 


0.32770666 . . . 


2.42990139... 


4.65118547... 






2*^ = 64 


0.32805801 . . . 


2.5381186... 


4.66503325 






2^ = 128 


0.32813334... 


2.4897987... 


4.66830065 . . . 






2^ = 256 


0.32814947... 


2.5099532... 
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0.3281538... 











